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Abstract 

Structural transitions in melts as well as their dynamics are con- 
sidered. It is supposed that liquid represents the solution of relatively 
stable solid-like locally favored structures (LFS) in the surrounding of 
disordered normal-liquid structures. Within the framework of this ap- 
proach the step changes of liquid Co viscosity are considered as liquid- 
liquid transitions. It is supposed that this sort of transitions repre- 
sents the cooperative medium-range bond ordering, and corresponds 
to the transition of the " Newtonian fluid" to the " structured fluid" . It 
is shown that relaxation processes with oscillating-like time behavior 
(u! ~ 10~ 2 s _1 ) of viscosity are possibly close to this point. 

1 INTRODUCTION 

It is known, that properties of both crystalline and amorphous substances 
often depend on their production conditions, and in many respects the ini- 
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tial liquid state determines these properties. At the same time, the liquid 
state of many substances is insufficiently investigated now. In recent years 
it has been found, that unusual properties and phenomena are character- 
istic for many fluids, but it is difficult to explain these properties within 
the framework of the elementary representations of the structureless liq- 
uids. One of the most interesting challenging phenomena being studied is 
the "liquid-liquid phase transition". At present the greatest progress is 
achieved in examination of network- forming liquids (such as water). Be- 
sides the papers, which inform about observation of similar phenomena in 
simple single-component systems, such as P [1] or C [2], have appeared. A 
possibility of existence of liquid-liquid transitions in the metal melts is also 
being discussed [3, 4]. 

Whereas at present the existence of liquid-liquid transitions in the glass- 
forming melts is beyond question, their existence in ordinary single-component 
metal liquids seems to be surprising. In many respects it is because of lack 
of direct structure observations. Nevertheless, there are already indirect 
evidences that point to the existence of appreciable changes in the melts 
structures. In particular, it is the appearance of "structural ffuid"-like be- 
haviour of some single-component melts at temperature decreasing, that 
appears as a little step change of liquid viscosity which goes on with the 
change of the viscosity activation energy [3, 5]. 

Besides, there are facts, which suggest that the process of establishing 
a thermodynamic equilibrium in the system close to this point can have 
a non-monotonic oscillating character [6]. This effect was observed on the 
polytherms of viscosity, surface tension and magnetic susceptibility of metal 
melts, and the low-frequency modes to ~ I0 -2 s~ x predominated in these 
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oscillations [6, 7] several times. For example, the time dependence of the 
kinetic viscosity of liquid Co at T = 1600 °C is presented in Fig. 1(b). These 
results were obtained in [8] by the torsional vibrations method in isothermal 
regime. As usual [8], the temperature was rapidly elevated and maintained 
at the given level. In the figure one can see considerable fluctuating of the 
measured quantities of viscosity. They have a stochastic character, and their 
amplitude exceeds the experimental error. It is important that the system 
temperature is close to the point of the step change of viscosity (Fig. 2). 
The dispersion of viscosity values in this point considerably exceeds the dis- 
persions of viscosity values at other temperatures. The usual variation of 
viscosity is Au ~ 0.7 • 10~ 7 m 2 /s for discussed experiments, but this disper- 
sion grows by 3-4 times near the vicinity of the temperature at which the 
viscosity polytherms jump are observed. The spectral-correlation analysis 
of these data was shown in the following papers [9, 10]. The results have 
confirmed the supposition that a pronounced low-frequency mode is present 
in these data (Fig. 3). 

Unfortunately, so far, the reliability of the above data has been ques- 
tionable and true conclusions about the character and mechanism of the 
observed relaxation processes have remained unsolved. It is connected with 
the difficulties of developing specific precise experimental techniques and 
processing the obtained data as well as with the lack of the appropriate the- 
oretical models. However, it is natural to suppose that these effects can be 
generated by some local structure modifications. Since the oscillation relax- 
ation processes are observed in the neighborhood of the step change of liquid 
viscosity we believe these effects to be connected with each other. Besides 
these vibrations look like as a fluctuations reinforcement near the phase tran- 
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sition point (Fig. 2). Therefore in this paper we assume that liquid-liquid 
transition point is root of both the step change of liquid viscosity and the 
oscillating relaxation. 

2 THEORETICAL MODEL 

It is known that thermodynamic systems are metastable at the temperatures 
close to critical ones. This metastability can be described by addition of the 
non-linear part to non-equilibrium thermodynamic potentials of the system. 
On the other hand, such systems undergo strong thermodynamic fluctua- 
tions as well as the effect of external "noise" because of both the influence 
of the experimental plant and thermal fluctuations. Thus, one can suppose 
that the combination of strong none-linearity of the metastable system and 
fluctuations is the reason of the oscillating relaxation, but the nature of 
metastability is clear only in the case of melts, whereas the nature of the 
liquid-liquid transition in liquid metals is still unclear. Therefore, we think, 
that the problem of theoretical description of the oscillating relaxation is 
reduced to two problems: first of all, it is necessary to define the nature 
of the liquid-liquid transition in the metal liquid and give its theoretical 
model. After that one can discuss the dynamic peculiarity of this model, 
and the chance of appearance of the low frequency mode in the oscillation 
spectrum. 

In most of theoretical papers liquid is treated as the one consisting of 
adjoining to each other and interacting with each other elementary local 
volumes which include atoms of several coordination spheres [11, 12, 13, 
14]. It is expected that each of these local volumes has only a few energy- 



4 



profitable configurations of the nearest order (local states). The long-range 
order is absent here as the local groups of atoms can be oriented differently 
with respect to each other. The topological structures of these locally- 
ordered formations have some symmetry. One can suppose that this local 
symmetry corresponds either to the crystal symmetry (solid-like), e.g. f.c.c. 
or b.c.c, or to the icosahedral symmetry (normal- liquid) . We believe that as 
the temperature decreases both normal-liquid and solid-like locally favored 
structures (LFS) [11] can form clusters. These clusters have a finite size and 
continuously transform, since the thermal fluctuations of the system lead to 
continuous destroying old bonds and to simultaneous forming new ones. 

Let us consider the case when the temperature of the melt is close to that 
of the liquid-liquid transition. The foolproof theory of the structural tran- 
sitions in liquids is not available for the time being. Therefore to describe 
the metastable properties of the close to this transition system we will con- 
sider the simplest model. As well as it was in [11] we will describe the close 
to liquid-liquid transition system as the system that is near the gas-liquid 
transition. We believe that the order parameter is connected with density. 
Therefore, using the density as the order parameter we introduce the follow- 
ing free energy, which governs ip fluctuations near the gas-liquid-like critical 
point T cr [15]: 



7(T, P) is the function of temperature T and pressure P, the parameters g 
and a are weakly depending on temperature. 
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3 DESCRIPTION OF THE RELAXATION 



Let us consider the fluctuations of the order parameter of metastable liq- 
uid in the region of phases co-existence (7 ~ 0). In order to describe the 
relaxation dynamics of this system we will consider the H-model [19]. This 
model is represented by the following equations system: 
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where r = a(T — T cr ), h and A are infinitesimal applied fields. The first 
equation describes the dynamics of the order parameter ip, and the second 
one describes the dynamics of the transverse part of the momentum density 
v, (P ± is a projection operator which selects the transverse part of the 
vector in brackets, A is the coefficient of self-diffusion, 770 is the viscosity, 
go = 1/ A770 is the mode-coupling vertex, the functions £ and 6 are the 
Gaussian white noise source: 

(6(x, t)8(x', t')) = -XV 2 5(x - x')8(t - t'), 

<&(x, t)Sj(xf, 0) = -??oV 2 5(x - x')S(t - t')5ij. 

The critical properties of this model are known [19]. In particular, it is well 
known that the viscosity of such a system depends on the character time 
of the experiment t c h = l/^ch- In the case of the results discussed above, 
this quantity is the frequency of torsional vibrations. Therefore, one can 



really expect that the change of the viscosity around T cr will be prolonged 
in some temperature interval T\ > T > T*(uj) rather than jump-like (T\ is 
the temperature of formation of the metastable phase). It agrees with the 
experimental observations. Note, that this interval depends on the frequency 
of torsional vibrations uj: 



{y and z are corresponding static and dynamic critical exponents). At the 
low- frequency uj < \jt c i (t c i is the life-time of the clusters ) the value of 
viscosity r)(u>) will be slightly dependent on the correlation length £ c (clusters 
size). But at the higher frequency uj » l/t c i ~ the response to the 
mechanical perturbation will be determined by the scale which is smaller 
than the cluster size, and it will not depend on the frequency. 



4 THEORETICAL EXPLANATION OF THE VIS- 
COSITY TIME OSCILLATING 



Usually model (2) is used for theoretical description of critical dynamics of 
gas-liquid transitions and transitions in binary fluids. The system inertia 
is not taken into account. The point is, the critical dynamics is usually 
investigated, and in this case the rescaling operator increases the importance 
of uj, relative to puj 2 , and after the renormalization procedure the term puj 2 
may be neglected. However, we will investigate the nonlinear stochastic 
system which is not found correct at the critical point. In this case the 
non-linearity of the system can lead to the state when even its weak inertia 
will essentially influence long-time dynamics [20]. In order to take it into 
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account, it is necessary to add the proportional to double t-derivation term 
to the first equation of the system: 



To analyze this stochastic model, one can employ the standard method of 
the stochastic deriving functional [21] and theory of perturbation. According 
to these methods, the correspondent field model will be described by a set 
of basic {ip, v} and supplementary {ip', v'} fields, and the effective action 
will have the form (Fig.4) of: 

5($) = -\ip'd 2 i>' + 4>'[- P d 2 4>-d t ^- Ad 2 (<9 2 V> - rV> - 7V> 2 - # 3 )- 

- vdip] + A- 1 ^ V0V + v'[-d t v + \' 1 g 1 d 2 v + 4>d(d 2 4>)]. 
The propagators of fields ip and v have the form of: 

<w> = <^v> T = — , m - 2U ' 2 lpl 



—uj 2 — iauj + Ek ' {—uj 2 — iaoj + £&) (— uj 2 + iau; + £&) ' 



(ot ) = (v v) = — — -J , (w) = 



-iwA# + fc 2 ' (-^A# + & 2 ) 2 ' 

where 

e fc = -k 2 (k 2 +r ), a = -, 
P P 

and ro is a renormalized quantity r, /c 2 is the impulse k squared. Below we 

will consider only the propagators of the conservative order parameter field 

ip. For calculation it is convenient to use the (k, ^-representation, in this 
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case, the propagators become: 
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when 4efc > a 2 . It is important that the effective viscosity depends on ui, 
and it can be represented as 



r/(r, = r/o 



where is a coupled-mode contribution to the response function, p is 
an external impulse. In an one-loop approach it can be represented in the 
diagram form (Fig. 5), in mathematic (k, ^-representation this contribution 
has the form of: 



T, v > v (t, p) 
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in case of 4e^ > a 2 . If we switch from integration at k to integration at 
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oj = y Ask — a 2 in limit t — ► oo and p — ► we can obtain: 

Jim p) = (pA) ^ f dw [frioj) sm(ujt) + / 2 (w)(cos(wt) + 1)] , 

p^o p* 16V2a 7 

where 
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Fig.6 shows the qualitative aspect of this spectrum. One can see tow peaks 
in this figure. The uj = peak corresponds to the Goldstone mode and it is 
explained by the presence of preservation law for ip and v. Another peak 



1 / 68.4Ap 

OJ — -4 / 1 



corresponds to the low-frequency oscillations of the system. Thus, it is pos- 
sible to conclude that the observed oscillations are noise-induced, analogous 
to the noise-induced oscillations in bistable systems, and excitation of the 
low-frequency modes in the system can be a noise-induced transition [22] 
(non-equilibrium phase transition). 
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Then the condition of observation of these oscillations is £ c ~ [68.4Ap] 1 , 
and one can estimate the correlation length quantity: Since A = D a a,Q, where 
D a ~ 10~ 9 m? I s is the diffusion coefficient, ao ~ 1CT 9 m is the interatomic 
distance, and p is the relaxation time, we believe that these values are of the 
order of the torsional vibration period (it is p ~ 1 s in our experiments [8]). 
Then the correlation length is £ c ~ 10~ 7 m. One can anticipate that it is the 
size of the flickering clusters of the metastable phase. It should be noted that 
this value is close to the size of the "Fischer clusters" which was observed 
in supercooled liquids [3] , but direct observation of these formations at the 
relatively high temperatures is not available for the time being. 

5 CONCLUSIONS 

It is well known that liquid is a nonuniform system, and its structure 
is characterized by the presence of the structure clusters in it. Recently 
the polymerization-like processes have been observed in metal systems by 
small-angle neutrons dispersion experiments [3]. In these experiments the 
snowflake-like large-size heterogeneity was discovered. We believe that sim- 
ilar processes of liquid-liquid transition lead to appearance of structured 
fluid (rheology) properties in liquids at relatively low temperatures, and to 
the jumps in the viscosity polytherms of some liquid metals. 

In an attempt to explain the observed oscillations of viscosity, we assume 
that presence of liquid-liquid transition point is the possible root of appear 
of the properties oscillating in the explored systems. We believe that they 
are caused by the exterior noise. Nonlinearity of the system in the region of 
its metastability is the reason of the increase of the fluctuations dispersion, 
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and the possible reason of the oscillating character of these fluctuations. The 
latter is caused by the low-frequency peak in the spectrum, and we believe 
that this effect analogous to well-known noise-induced transitions in bistable 
systems [22]. We would like to mark that one ought not to consider the 
period of the viscosity oscillation as a life time of the metastable subsystems 
(clusters). The sizes of the such subsystems are about correlation radius 
~ £ c and their fluctuations frequency is a reciprocal value for the lifetime 
tcl ~ £c ~ 10~ 10 s. As indicated above the viscosity oscillation is the 
dynamic effect, which is determined by influence of the fluctuating moving 
to the hydrodynamic moving, and inheres in the whole system. 

Further examination of nontrivial dynamic properties of such systems 
will allow us to confirm or discard our hypothesis. 

This study was supported by the RFBR grant (01-02-96455 r2001ural). 
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Figure 1: Time dependences of kinematic viscosity of liquid copper [9, 10] (a) 
and liquid cobalt (1870 K) [8] (b) close to the temperature of step change 
of viscosity. 
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Figure 2: Temperature dependence of kinematic viscosity of liquid cobalt 

(a) and time dependence's of liquid cobalt viscosity at various temperatures 

(b) . 
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Figure 3: Time dependences of autocorrelation function for the liquid copper 
(a), and liquid cobalt (b) viscosity close to corresponding temperatures of 
step change of viscosity. 
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Figure 4: Graphs a,b,c,d correspond to propagators < ipip >, < ipip' >, 
< vv >, and < vv' > accordingly, graphs f,g,h correspond to vertexes 
v'ipd(d 2 tp), Xgtp' d 2 (ip 3 ) , and ip'vdip accordingly. 
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Figure 5: The graph corresponds to the response function E</„ in the frame- 
work one-loop approximation. 
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Figure 6: The spectrum of viscosity oscillations (f(u>) = y / ' f\(u) 2 + /2(w) 2 ) 
in linear-linear (a) and logarithmic-linear (b) coordinates. 
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